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CALCULATED EFFECTS OF THE LATERAL ACCELERATION 
DERIVATIVES ON THE DYNAMIC LATERAL STABILITY 
OF A DELTA-WING AIRPLANE 
By John P. Campbell and Carroll H. Woodllng 

SUMMARY 


Calculations have been made of the dynamic lateral stability of 
a 6o° delta-wing interceptor airplane with the lateral acceleration 


derivatives 


s 


and C^^ included and neglected. Calculations were 


made for angles of attack of 10°, 20°, and 50 °, with the airplane flying 
at sea level and at an altitude of 50,000 feet. Including the lateral 
acceleration derivatives in the calcvilations caused changes in stability 
that were small at 10° angle of attack where the values of these deriva- 
tives were small, fairly large at 20° angle of attack where the deriva- 
tives were larger, and very lfiu:ge at 50° angle of attack where the deriv- 
atives were very large. These results Indicate the necessity for including 
the lateral acceleration derivatives in calculations of dynamic lateral 
stability. In practically all cases, including these derivatives caused 
increases in the damping of the Dutch roll oscillation. The effects 
of . Cjip 8Ji<i 0^^ varied greatly when the airplane was assumed to have 

different values of static directional stability 


and effective 


dihedral effects on stability of varying th4 yawing deriva- 

we*re different from the effects of varying 0^* 


tives C 


n 


and C^' 


and Cl 




r *'r 

especially at the high angles of attack. 


INTRODUCTION 


Recent NACA research with wind-tunnel oscillation testing equipment 
has shown that large values of the lateral acceleration derivatives Cn" 


3 


and C^A sometimes are obtained with highly swept vm: 
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of attack. (See refs. 1 and 2.) In the past, these derivatives have 
usually been neglected in making dynamic-lateral-stability calculations 
for airplanes because little information was available to permit rea* 
sonably accurate estimates of the derivatives. Now that experimental 
data on Cn* an<i have become available, a theoretical investiga- 


tion of the effects of these derivatives on dynamic lateral stability 
has been \ondertaken. In this investigation, calciaations are being made 
for a variety of configurations and flight conditions with the lateral 
acceleration derivatives .included and neglected. Some preliminary results 
of this investigation, consisting of stability calculations for a delta- 
wing interceptor airplane, are presented in this report. 


Since the values of Cn^ and appear to be appreciable only 

at moderate and high angles of attack, the calculations for the delta- 
wing interceptor were made only for angles of attack of 10*^, 20®, and 50 • 
The calculations were made both for sea level and for an altitude of 
50,000 feet." The airplane was assxmied first to have values of Cn^ 


and 


of zero and then to have values of these derivatives similar 


to those obtained from oscillation tests of a delta-wing model in the 
Langley free -flight tunnel. 


SYMBOLS AND COEFFICIENTS 


(f) angle of bank, radians 

i|f angle of azimuth, radians 

3 angle of sideslip, v/V, radians 

V airspeed, ft /sec 

V sideslip velocity along Y-axis, ft/sec 

V sideslip acceleration along Y-axis, ft/sec^ 

p rolling velocity, d^/dt, radians/sec 

r yawing velocity, dijr/dt, radians/sec 

3 rate of change of angle of sideslip, v/V, radians /sec 

oib/2V reduced-frequency parameter 
(D ang ular velocity, radians/sec 
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a 


3 


P 

q 

b 

s 

w 

m 

g 


€ 


a 


7 






Kz 


mass density of air, sl\igs/cu ft 
dynamic pressure, ^pV^, Ib/sq ft 

wing span, ft 

wing area, sq ft 

weight of airplane, lb 

mass of airpleine, w/g, slxags 

acceleration due to gravity, ft/sec^ 

relative-density factor, m/pSb 

angle of attack of principal longitudinal axis of inertia, deg 
(see fig. 1) 

angle between reference axis and principal longitudinal axis of 
inertia, deg (see fig. 1) 

angle of attack of reference sixis, deg (see fig. 1) 
angle of climb, deg (see fig. 1) 

radius of gyration in roll about principal longitudinal axis of 
inertia, ft 

radius of gyration in yaw about principal normal axis of inertia, 
ft 

nondimens ional radius of gyration in roll about principal longi- 
tudinal axis, 

nondimensional radius of gyration in yaw about principal vertical 
axis, kzo^b 

nondimensional radius of gyration in roll about longitiodinal 
stability axis, 

nondimensional radius of gyration in yaw about vertical stability 
axis, \ /k 2 ^^cOs^i 1 + Kxp^sin^T] 


n/ 


KZo^sin2T) 
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go: 



nondimens ional product-of -inertia factor, 

^Note that negative for positive values of t).^ 

trim lift coefficient, Lift/qS 

rolling-moment coefficient. Rolling moment/qSb 

yawing -moment coefficient. Yawing moment/qSb 

lateral-force coefficient. Lateral force/qS 

effective -dihedral derivative, rate of change of rolling -moment 
coefficient with angle of sideslip, dG^^^P, per radian 

static directional- stability derivative, rate of change of 
yawing-moment coefficient with angle of sideslip, ^Gn/^P, 
per radian 

rate of chsmge of lateral-force coefficient, with angle of side- 
slip, ^Gy/^P, per radian 

rate of change of rolling-moment coefficient with lateral- 

acceleration factor, per radian 

I 2V 

rate of change of yawing-moment coefficient with lateral- 

acceleration factor, ^Gn/^^. per radian 

I 2V 

rate of change of lateral-force coefficient with lateral- 

acceleration factor, per radian 

2V 


damping-in-yaw derivative, rate of change of yawing-moment coef- 
ficient with yawing-angular -velocity factor, per 

7 2V 

radian 

rate of change of yawing-moment coefficient with rolling-angtilar- 
velocity factor, per radian 
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Cl. 


Cl, 


CYr 

t 

s 

D 

P 

'^1/2 


' 0 

V i 


rate of change of rolling-moment coefficient with rolling- 

angular -velocity factor, per radian 

7 2V 

rate of change of rolling-moment coefficient with yawing-angular 

velocity factor, per radian 

7 2V 

rate of change of lateral -force coefficient with rolling -angular 

velocity factor, bCyl^^, per radian 

7 2V 

rate of change of lateral-force coefficient with yawing -angular- 
velocity factor, per radian 

time, sec 

nondimens ional time parameter based on span, Vt/b 
differential operator, d/ds 
period of oscillation, sec 

time for amplitude of oscillation to change by a factor of 2 
(positive value indicates a decrease to half -amplitude , nega-- 
tive value indicates an increase to double amplitude) 


A,B,C,D,E coefficients of lateral-stability equation 


EQUATIONS OF MOTION 


The nondimens ional linearized lateral equations of motion, referred 
to the stability axes, used to calculate the stability roots are as 
follows; (These equations are the same as those of reference 3 except 
for the axMition of the ^Ip^ ^Yp ieritis and a change in the 

sign of Kxz') 


Rolling moment 

2^x(Kx2D2^^ - 


= Cip3 + I C2*D3 + I i Ci^Di/ 
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Yawing moment 

- %zD^) » CnpS + I Cn-DS + | 0„pll{l + | Cn^D* 

Lateral force 

2 ^^D 3 + D\trj = CYp 3 + ^ Cy^D 3 + i ■*■ ■*■ 2 ■*■ ^^’1' 

When 0 q 6 ^® is substituted for (^, i|foe^® for \|r, 3 q6^s for 3 in 

the equations written in determinant form, A mvist be a root of the 
stability equation 

AA^ + BA^ + CA^ + DA + E = 0 

where 

A = 8h5^Kx%z^ - Kxz^) - 2 h2(kx%z^ - Kxz^)cy^ 

B = -2h(2mKx%22cyP + mKx^C^^ + “ SMKxz^Cy^ + dKxzCj^ + 

MKxzCnp - ^ Kx^Cn^Cy^ - mKx^Cq^ - i Kz^Ci^Cy- - i KyzCl^Cy^ - 
^J^^XZC^^ - ^ Cn^KyzCy^ + ^ Kyz^Yp^n^ + J ^z^YpCl^ + 


i Kx^Cy^Cn- + i KxzCy^C^^ 
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C = MKjj^nr^Yp + + MK^^C^pCyp + | + M^zClr^Yp + 

l^^l%xzC^p + ^^Cnp%Z^3 “ 2 ~ " 

MKz^Yp^Zp ■ ^^X^CY^<^np " MKxzCY^^lp “ | ^n^Plp^Y'^ “ 


I nCipCn- + I Cn^Ci^CY- + ^ nCn^Ci^ - nCj^CxzCnp - 

7 - uKxzCzpCL tan 7 + 


i C, Cv C„; - i C„ Cv C 7 I - i C 7 Cv CnA + - Cn Cy Ct- 
3 Ip Yr H 3 3 np Yr ip 3 ir Yp np 3 % *p ^p 


/^ 


D = - i Cnj.CipCY 3 - ^CZpCnp + i CnpCz^CYp + ^CnpCZp - 


2^lC];KxzCnp " tan 7 - 


^MKxz^ZpCL tan 7 + i C2^CnpCY^ - ^ ^n^C^pCY^ - 


— C7 Cj» Cy + — C-, C7 Cy + — Ct C,~ Cl* ~ — CtCi Cj-* + 
ir “p ^p 1}. ^r ^p ^p 4 % ^p 1). “3 


i CtCt C_. tan 7 - ^ CrC„ C7. tan 7 

l^. ^ ip np ’ li ^ np ip 


E = I Cl(c„^C ,|3 - Ci^C^p) + i Cl tan r(cipC„j, - CnpCj^^ 
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CALCULATIONS 


Calciilations were made to determine the period and damping of the 
oscillatory mode and the damping of the aperiodic modes using the eqvia- 
tions presented in the preceding section. The dimensional and mass char- 
acteristics assumed for the airplane are presented in table I and the 
flight conditions for which calculations were made are given in tables II 
and III. Calculations were made for angles of attack of 10 °, 20 °, and 50 ° 
for sea level and for an altitude of 50,000 feet. 

Values of the stability derivatives used in the calculations are 
given in tables II and III and plots of variations with angles of attack 
of some of the derivatives are shown in figure 2 . For 30 ° angle of attack, 
four different combinations of the important sideslip stability deriva- 
tives Cn and Cj were assumed as follows: 


Combination 

Cn 3 

r 

A 

-0.0573 

0 

B 

-.0573 

-.0573 

C 

.0573 

0 

D 

. .0573 

-.0573 


Studies of force-test resxilts for several delta-wing configurations have 
indicated that, depending upon the particxilar configuration, any one of 
these combinations of Cn„ and C 2 might exist at 30° angle of attack. 

p 3 

Values of the rolling derivatives C^p and Cn^ and the yawing deriva- 
tives Crip and Cj^ were estimated from experimental data obtained by 
the NACA on delta-wing configurations; and the derivatives Cyp^ Cy^, 


and Cyp 


were assumed to be zero. 


Values of and Cjp were estimated from the data presented 

in references 1 , 2 , and 4 . In maiking the estimates of C^^> values 
of Cn obtained from curved-flow tests (ref. 4 ) were subtracted from 
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the veilues of 
of Cn^ - Cn^; 


Ciir “ Cnp 
Clp; 


of reference 2. For estimating veilues 

and Cj measvired about the stability axes 


and values of Cn^, 



cos 0 


and 


+ C?* sin 0 
P 


measured about the 


body axes (ref. 1) were substituted in the equations relating the stability- 

axis and body-axis damping derivatives. The equations were then solved 

for C? • . The data of references 1 and 2 were obtained at values of the 

ip 


reduced-frequency parameter oi)/2V of about 0.21 to O.25. Subsequent 
oscillation tests of a delta-wing model (resiilts unpublished) have indi- 
cated large effects of frequency, particularly at angles of attack above 
about 15° ♦ The effects of frequency are such that greater values of Cn^ 

and Cjp are obtained at the smaller values of ai>/2V. In addition to 


these results, some results (also unpublished) have recently been obtained 
in the Langley stability tunnel which show that the values of Cn^. and 

Ci^ for a 60° delta wing in an oscillation are actually much greater 

them the values obtained in the curved-flow tests of reference 4 . The 
values of Cn^ and used in the present calculations are therefore 

probably smaller than they should be. These data also indicate that the 
values of Cn^ and used in the present calculations are too large 


since they were obtained by siibtracting the curved-flow-test values 
of Cn and C7 from the oscillation-test values of Cn„ - CnA 


and C^j. - Cjp. The overall significance of these additional unpublished 

results in connection with the present calculations will be discussed in 
the "Results and Discussion" section. 


Calcinations were made for each condition with Cn- and Ci* 

Up tp 

neglected, with Cn* included, and with both Cn^ and Ci^ included. 

In addition to these three calculations for each condition, two other 
calcvilations were made with Cn^ and Cjp neglected - one in which Cn^, 

was increased by an amount which would correspond to Cn^ - Cn^ and the 

other in which Cn^. and Cj^. were increased to correspond to Cn^ - 

and Cjj, - Cj^, respectively. These latter calculations were made to 

determine whether- Cn^. and C^^ might have the same general effects on 

stability as Cnp and ^2^- It might be reasoned that such could be the 
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case, since it appears that most of the terms of primary importance con- 
taining these four derivatives in the coefficients of the lateral-stability 
equation can be combined into terms involving - Cnp and 


RESULTS AND DISCUSSION 


The resiilts of the calculations are presented in tables II and III. 

Some of the results showing the effects of Cni and Ci* on the period 

P *'p 

and damping of the Dutch roll oscillation are presented in figures 3 
and 4. 


Effects of Cnp and 

The data of table II and figvire 3(a) show that for 10° angle of 
attack the effects of Cn^ and C^p on either the Dutch roll oscillatory 

mode or the aperiodic modes were quite small for both the sea-level and 
50,000-foot-altitude cases. The Dutch roll damping was increased slightly 
by including Cnp and still fvirther by including both Cnp and 

The relatively small effect of the lateral acceleration derivatives at 
10° angle of attack can be attributed to the fact that these derivatives 
are small at this angle of attack as shown by figxire 2. 

For 20° angle of attack, the data of table II and figure 3(b) show 
that substantial improvements in Dutch roll damping at sea level and at 
50,000 feet were obtained by including Cnp and Cjp in the calcula- 
tions. The changes in the dajnping of the aperiodic modes resulting from 
inclusion of Cnp and C^p at this angle of attack, however, were very 

small. 

The data of table III and figure 4 show that at 30° angle of attack 
the effects of Cn^ and C^^ were very pronounced for both the Dutch 

roll and aperiodic modes of some of the combinations of C^p and C^^. 

In most cases the damping of all the modes was increased. An extreme 
change from a high degree of instability to a high degree of stability 
of the Dutch roll oscillation was obtained with combination D by 
including the acceleration derivatives. On the other hand, including 
Cnp and C^^ in the case of combination A caused only relatively minor 

changes in Dutch roll damping. These results indicate that the effects 
of Cnp and Cj^ are greatly dependent on the values of the other 
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stability parameters such as Cnp and A detailed analysis of the 

results of this and other similar investigations might reveal the reasons 
for some of the differences in the effects of Cn^ and for differ- 

ent cases but more than likely it will not be possible to make any gen- 
eralizations regarding the effects of these derivatives. Previous studies 
of dynamic lateral stability have indicated that it is unwise to attempt 
to generalize regarding the effect of any given stability derivative 
because of the strong interdependence that exists between most of the 
derivatives . 


As pointed out in the Calciilations section, values of 
and Cjp used in the present study were estimated from oscillation test 


data obtained at values of oi)/2V from about 0.21 to 0.25* The results 
of the calculations, however, show that for the flight conditions con- 
sidered in this investigation the values of oib/2V are in most cases 
much smaller than 0.21. As stated previously, unpublished test results 
obtained with a delta-wing model subsequent to the tests of references 1 
and 2, have indicated that the values of Cn^ and Ci^ are much larger 

at the lower values of aib/2V, particularly for angles of attack above 15° 
On this basis it would appear that the values of Cq* and used in 


the present calc\ilations for 20° and 50° angle of attack are too small. 
On the other hand; as pointed out in the "Calculations" section, the 

and Cl indi- 


additional impublished resvilts on the derivatives 


^nj. 


cate that the values of Cnp and C^* \ised in the calculations are too 


3 


large. A preliminary analysis of the limited amount of data available 
at present indicates that these factors are roughly compensating so that 
the values of Cn^ and Cj^ used are at least of the right order of 


magnitude . 


In any event, it is unlikely that the changes in 



and Cjp which might be involved woxild alter the principal conclusions 

drawn from these calculations - that the effects of Cn^ and Cj^ are 

appreciable and that these derivatives should be considered in studies 
of dynamic lateral stability. 


Effects of C. 


nv 


and C-j 


The results of calculations made to determine whether Cn and C? 

•‘‘r ‘'2 

have the same effects on stability as Cn^ and Ci^ are presented in 
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tables II and III. In each group of five rows in 
second and fourth rows afford a comparison of the 


these tables, the 
effects of Cn* 


and Cv, while the third and fifth rows give the results when either 
nr 

Cnp aJ3d Cjp or Cn^ and were changed. 


It is apparent from the results presented in tables II and III that 
the effects of Cn^ and cannot always be properly simulated by 

making corresponding changes in Cn^. and The differences between 

the effects of the p and r derivatives were especially great at 30° 
angle of attack. Since, for the purpose of this report, this general 
result is the only one of interest in connection with the and 

calculations, no detailed discussion of these results will be given. 


CONCLUDING REMARKS 


The results of the dynaraic-lateral-stabillty calculations for the 
6o° delta-wing interceptor airplane with the lateral acceleration deriva- 
tives Cn^ and C^^ Included and neglected can be summarized as follows 


1. Including C^^ and Cj^ 


in the calculations caused changes in 


stability that were small at 10° angle of attack where the values of 
these derivatives were small, fairly large at 20° angle of attack where 
the derivatives were larger, and very large at 30 ° angle of attack where 
the derivatives were very large. These results indicate the necessity 
for including the lateral acceleration derivatives in calculations of 
dynamic lateral stability. 


2. In practically all cases, including Cn^ and caused 

increases in the damping of the Dutch roll oscillation. 


3. The effects of and varied greatly when the airplane 


was assumed to have different values of Cn^ and 
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4. The effects on stability of varying Cn^. and Ci^ were differ- 
ent from the effects of varying Cn^ and especially at the high 

angles of attack. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , November l6, 1954. 
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TABLE I 


DIMENSIONAL AND MASS CHARACTERISTICS ASSUMED FOR AIRPLANE 


Weight, lb . . . . . . 
Wing loading, Ib/sq ft 

(l^)h=0 


K. 


■Xo 

KZo^ 


€, deg 


22,850 

34.5 

11.85 

0.0155 

0.0844 


1.2 


Wing : 

Area, sq ft 662 

Span, ft 58-1 

Aspect ratio 2.2 

Sweepback of leading edge, deg 60 

Vertical tail: 

Area, sq ft 6.8 

Span, ft 8.7 

Aspect ratio 1»1 

Sweepback of leading edge, deg 60 





CALCULATION CONDITIONS AND RESULTS FOR 50° ANGLE OF ATTACK 
















TABLE III.- Concluded 
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Cnp Cip 

A -D573 D 
B -.D573 rD573 
C D573 D 
D D573 -.0573 


Tl 1 00 
2 


Periods >20 seconds are indicated 
by numbers beside symbols. 

A long period oscillation was also 
obtained with combination C.{SeetableIlD 


12 

P. sec 


20 > 20 


(a) Sea level (h = O). 

Figure 4.- Effect of Cn^ and C^* on period and dan5)ing at 50° angle 
of attack for various combinations of Cno and C7 . 
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